Radiation reaction for multipole moments 
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We propose a Poincare-invariant description for the effective dynamics of systems of charged 
particles by means of intrinsic multipole moments. To achieve this goal we study the effective 
dynamics of such systems within two frameworks - the particle itself and hydrodynamical one. 
We give a relativistic-invariant definition for the intrinsic multipole moments both pointlike and 
extended relativistic objects. Within the hydrodynamical framework we suggest a covariant action 
functional for a perfect fluid with pressure. In the case of a relativistic charged dust we prove the 
equivalence of the particle approach to the hydrodynamical one to the problem of radiation reaction 
for multipoles. As the particular example of a general procedure we obtain the effective model for 
a neutral system of charged particles with dipole moment. 

PACS numbers: 03.50.De, 47.75. +f, 52.27.Ny, 04.20.Fy 



I. INTRODUCTION 

We address this paper to the problem of the construc- 
tion in the classical electrodynamics framework of an ef- 
fective model describing in some approximation the dy- 
namics of an extended object consisting of charged parti- 
cles. The word effective means that we take into account 
the self-interaction or radiation reaction force exerting 
on such a system and thereby eliminate from dynamical 
equations the electromagnetic field produced by charged 
particles. 

The simplest models of this kind are the model of a 
point charged particle obeying the Lorentz-Dirac equa- 
tions m, y and its generalizations to curved or higher 
dimensional space-times [1,11,1^. Other models, which 
can be regarded as ad hoc effective models of this kind, 
emerge at the classical description of charged spinning 
particles H 0, B H E lli ■ In [13 the effective model 
is formulated and investigated for a high-current beam 
of charged particles taking into account the leading self- 
interaction correction. In the paper [l3j the effective 
equations of motion are obtained for the center of mass of 
a rigid charged body of arbitrary shape at the pointlikc 
limit. As expected they are the Lorentz-Dirac equations. 
A nonrelativistic approach (1/c expansion) to the effec- 
tive dynamics of a system of charged particles considering 
the radiation reaction due to electric (magnetic) dipole 
and quadrupole radiation can be found in 14, 15]. 

In the present paper we obtain a Poincare-invariant ef- 
fective model for a system of charged point particles as 
the model of a point particle with some internal degrees 
of freedom like the intrinsic dipole moment, the intrin- 
sic magnetic moment etc. Briefly, the procedure is as 
follows. We solve the Maxwell equations for arbitrary 
worldlines of constituent charged particles. Then we find 
the mean electromagnetic field in a small neighborhood 
of the system under consideration by averaging over all 



its small-scale fluctuations^ and substitute this field into 
the expression for the Lorentz force. Using the obtained 
system of equations we derive the equations of motion 
for the center of mass of the object and the evolution of 
intrinsic multipole moments. Neglecting terms of higher 
order of smallness we break the chain of equations and 
arrive at the closed system of equations of motion for the 
effective model. These steps are scrutinized in Sections 
lillandlml 

In Section ^ we also give a Poincare-invariant defi- 
nition of the intrinsic multipole moments. Besides, as 
the particular case of a general construction the effective 
model for a neutral pointlike system of charged particles 
is obtained at the end of Section IIIII where we estimate 
the energy losses of such a system and investigate its free 
dynamics. 

In Section IIVI we consider an alternative approach to 
the description of systems of charged particles. Namely, 
we pass from the particle framework to hydrodynamical 
one and suggest an obvious generalization of the stan- 
dard action functional for particles to the hydrodynam- 
ical case. This action is found to describe a relativis- 
tic dust. Then we generalize the definition of intrinsic 
multipole moments to the hydrodynamical approach and 
prove the equivalence of particle and hydrodynamical ap- 
proaches to the problem of radiation reaction for multi- 
pole moments. 

In addition, in Section IIVI we generalize the action 
functional for a relativistic dust to a relativistic per- 
fect fluid with pressure and define in a Poincare-invariant 
manner the intrinsic multipole moments for extended ob- 
jects (charged fluids or systems of charged particles) ap- 
proximated by branes, i.e. we give a relativistic deflnition 
of the linear density of dipole moment for a string or the 
surface density of quadrupole moment for a membrane 
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^ The averaging procedure is realized as the regularization of self- 
force and analogous to averaging procedures applicable in t he 
mean field and renormalization group methods (see, e.g., Il6l ). 



etc. In concluding section we summarize the main re- 
sults and outline the prospects for further investigations. 



II. EQUATIONS OF MOTION AND 
MULTIPOLES 



In this section we recall some basic formulas concerning 
classical electrodynamics of many-particle systems and 
define multipoles for such systems in a Poincare- invariant 
manner. 

Let R'^'^ be 4-dimensional Minkowski space with co- 
ordinates {a;'^}, n = 0, 1,2,3, and signature (-I-, — , — , — ). 
In the space-time given a system of N electrically charged 
point particles with trajectories Xa{Ta), a = 1,...,N. 
The dynamics of the system in question are governed by 
the action functional^ 

S[Xa{Ta),A{x)] = ^ -TOa / dTa\/il - / d'^xA^.j^ 



a=l 

1 
16^ 



(X X-tniy-t , 



(1) 

where Fp,y = 9[^A,^] is the strength tensor of the electro- 
magnetic field, overdots denote derivatives with respect 
to the parameter on the worldline and the electric cur- 
rents of point particles are 



Jaix) = ea dTa6'^{x - a;a(Ta))i^ (t^) . 



(2) 



The equations of motion for the action functional (^ 
constitute the system of Maxwell-Lorentz equations and 
in the Lorentz gauge d^A^^ = are given by 



d 

' dTa 



N 



= e,n(^aK, n^'^ = 4^^j^ 



(3) 

To obtain the effective model for a system of charged 
point particles we, first of all, should solve the Maxwell 
equations for arbitrary worldlines x^{Ta)- The casual so- 
lution to the Maxwell equations can be constructed by 
the use of retarded Green's function G{x) associated with 
the d'Alembert operator |17l |: 



G{x) 



27r 



6{x^ 



Thus the Lienard-Wiechert potentials 

N 
A^^{x)=2Y,ea dTa0{x° ~ X°{Ta)) 



a=l 

X S{{x - 



(4) 



(5) 



Xa{ra)r)x^a{ra) 



give the solution to the Maxwell equations. They are 
interpreted as the electromagnetic field created by the 
system of charged particles. A solution of the homo- 
geneous Maxwell equations, which can be added to the 
Lienard-Wiechert potentials, is regarded as an external 
field. 

For our purposes it is useful to parametrize the tra- 
jectories of particles in the following way: let z^(t) be 
a naturally parametrized worldline in Minkowski space, 
then we parametrize the trajectories x^{Ta) by the pa- 
rameter r and claim that 



ipMC'M = 0, i^{T) = 1, 



(6) 



where $^(t) = x^{t) — z^^{t). The gauge © properly 
fixes parametrizations on the worldlines under the con- 
dition 



Zpi'a > 0. 



(7) 



which is obviously satisfied. The worldline z^{t) is de- 
fined by the requirement 



N 

E 



rua^ilir) = 0, 



(8) 



and we call it as the center of mass. 

By the intrinsic electric and magnetic multipole mo- 
ments of the system of point charged particles we under- 
stand irreducible components of the tensors'^ 



AT 

E 

a=l 



'^aS/ji ■ ■ ■ S/i„ J 



N 

E 

o=l 



^aQ,---^i,.J>^: 



(9) 



respectively, where prj;^ 



S'. 



z-z^ 



Evidently, they are 



orthogonal to i'^. For example, the Lienard-Wiechert 
potentials © can be rewritten in terms of the multipole 
moments as 



A^{x) = qDz^ + D[d^- ripd^z^] 

+ ^D[QDz^ + XPQp^X^D^z^ + S^^n" 
+ izP-nP)Qp^{S:^n'^ip)] + 0{e), 



(10) 



where one must assign r = Tret after all differentiations. 
We have introduced the notations 



^ fi — -^/i 



D = 



Rdr' 



z^i-r), 



X' 



R = ipXP, 



= 0, 



Xp 
R ' 



(11) 



^0 !.=.„, > 0, 



•^ Hereinafter we use the system of units in which c = 1. The 
Greek indices are raised and lowered using the Minkowski metric 
•qfxv on M?'^. Square (round) brackets at pair indices denote 
antisymmetrization (symmetrization) without one half. 



^ Usually electric multipole moments are defined as traceless ten- 
sors, since at large distances in stationary case the electromag- 
netic field of a charged object depends only on their traceless 
parts. But in the dynamical case the whole tensors should be 
employed. 



and the multipole moments 
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/v 



'? = X1^'^' '^M^XI^'^C' 



w 



a=l 



(12) 



Qp.v — 2_^ ^ai^C,^, Sfii, — 2_^ ^ai[^ip Pr^], 



a=l 



a=l 



which are the total charge, intrinsic dipole, quadrupole 
and magnetic moments respectively. Besides Q is the 
trace of the quadrupole moment. Higher terms in the 
expansion (|10|l . which we denote as 0(C'^), can be ex- 
pressed in terms of the multipole moments © as well. 

In particular, when all the multipole moments are con- 
stant and the system moves freely we obtain a relativistic 
generalization of the well known expression |l4| : 






^m(^) — "d^P^' 



2i?3 



where Q^^ = Q^^ - \ pr^^ 



(13) 



III. EFFECTIVE DYNAMICS OF MULTIPOLE 
MOMENTS 

In this section we substitute the Lienard-Wiechert po- 
tentials into the Lorentz force and regularize the resulting 
ill-defined expression. After the regularization we obtain 
an infinite series in a regularization parameter and ^°. 
To truncate this series we impose constraints on charac- 
teristic scales of the charged object under consideration. 
Thereby we derive the equations of motion for the ef- 
fective model associated with a system of point charged 
particles. In particular, we obtain the effective model for 
a neutral pointlike object with intrinsic dipole moment, 
estimate the energy losses of such a system and describe 
its free dynamics. 

The effective equations of motion of charged particles 
at lower orders in ff, look like 



d_ 
'd7 






(14) 



where F^i, is the strength of the external electromag- 
netic field taken at the point z^(r). As one can see the 
linear order in ^ of particle's momentum is orthogonal to 
i^, that is the reason that we define the intrinsic mag- 
netic multipole moments as projected by pr^. The field 
strength F"[, is constructed from the potentials ((SJ and 
is given by 



N 



F;:{xa{T)) = 4^6,/ dse{x' + c°,) 

h=l -^ 



b=i - (15) 
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FIG. 1: The part of an infinite lattice depicting contributions 
of the integral 1151 . The vanishing terms are denoted by crosses. 
They are the terms at negative integer powers of e and the terms 
at g'^e 2 , k < I. The dotted triangle singles out the contribu- 
tions, which we take into account for a charged system meeting 
1191 . Three marked dots depict contributions of the self- force to 
the equations of motion for the center of mass of a neutral pointlike 
object to the accuracy of the first radiation correction. 



where X^ = z^(t) - z,{s) and C^' = ^^ir) - C^(s). Now 
we expand the integrand of H15|l in powers of S, and arrive 
at 

6'ix')iX[^ + 0ii^ + s"ix'){{e'r + 2^"^;') 

y<{^^;;+X[,)xl^+2d'"{x'){Cfx''f{^f^+Xi,)il^ + ... 

(16) 

Here dots denote neglible in our approximation terms 
(see below). 

After integration every term in the series (|16l) gives rise 
to infinities, because the ^-function and its derivatives 
are ill-defined at the vertex of light-cone (see, e.g., [ig 
and discussions in |a,llMl23l)- To cure this problem we 
have to regularize the integrals, i.e. to represent them 
as sequences of converging integrals. We apply the so- 
called "point-splitting" regularization of the (5-function 
and its derivatives which lies in substraction of a posi- 
tive number e (the regularization parameter) from their 
arguments. Under this procedure the support of the 6- 
function transforms into a hyperboloid which is a smooth 
manifold. That is why the integrals become converging. 

In fact, the regularization of the (5-function is the regu- 
larization of the Green function Q . Therefore the phys- 
ical meaning of the regularization procedure consists in 
averaging over all small-scale fluctuations of the electro- 
magnetic field produced by the system of charged parti- 
cles up to the scale e. That is to say the regularization 
parameter is a characteristic scale of fluctuations of the 
electromagnetic field of the charged object. 

A useful mathematical framework for handling inte- 
grals with integrands like ((TC|l is elaborated in |23 and 
we do not enlarge on it here. The structure of contribu- 
tions of the integral (|15|l is as shown in Fig. Q] All the 
terms at half- integer powers of e are Lagrangian [20|, i.e. 
can be obtained by varying an effective action, which in 



turn is obtained from the action functional Q by sub- 
stituting the Lienard-Wiechert potentials lO in it and 
applying the regularization procedure. The terms at in- 
teger powers of e are responsible for radiation losses and 
not restored by the effective action. 

Let I be a characteristic scale of variations of the trajec- 
tory z(t). Then, firstly, we average over* all oscillations 
of the fields ^'t(j) the frequencies of which are greater 
than l^^ . In that case the n-th derivative of ^"(r) with 
respect to r is of the order ^/Z", where ^ is a character- 
istic scale of fluctuations of the variables ^° (r) . As the 
result for the integral (|15l) we have an infinite series in 
two dimensionless variables ^jl and e^l"^ jl. 

Secondly, one can see from Fig. ^the necessary condi- 
tion for the obtained asymptotic expansion makes sense 
is 



i < e^/2 < /. 



(17) 



Physically, the condition (|17|l means that we consider a 
charged system which is much smaller than the charac- 
teristic scale of fluctuations of the electromagnetic field 
(after averaging) while this characteristic scale is much 
smaller than the characteristic scale of variations of the 
trajectory zir). 

Thirdly, in the lack of any a priori data on the system 

* It is necessary so as to the time-scale of variations of the fields 
5JJ(r) would be larger than e^'^ . Otherwise the terms at higher 



of charged particles in question we make the assumption 
that there exists n such that 



(n+fc) 

M 






(18) 



in) 



for any positive integer fc. Here M denotes a magni- 
tude of the n-th multipole moment ©. In other words 
we assume that beginning from the n-th order multipole 
collective phenomena such as spontaneous magnetization 
etc. are absent in the charged object, i.e. there is no any 
high ordered structure in it. If that is the case then the 
leading contribution to the self-force is contained in the 
first n orders in ^ of the asymptotic expansion. 

In subsequent calculations we mostly deal with the case 
n < 2, moreover we take 



e//«(ei/70^ e/f<^/l. 



(19) 



These conditions mean that we consider a charged object 
which is much smaller than e^", but is not so small to 
be pointlike. Some of the above requirements could be 
relaxed, however it deserves individual investigations. 

So, neglecting the terms of higher order of smallness 
we obtain 

powers of e make greater contributions than at lower ones. 



J 



F: 



N . 5 _3 

6=1 '^ 



16 



(2 - Sz^ep"" - Qi'Ozipi.] + 2z[^it] + 2eti.] - 2^[';z-,] - -^f^i.] 



£2 



6^[/^^-l+32 



(4) ... 5 2 

3 z [^z^] + 2 z [^Zj,] + -z Z[^z^^ 



, (20) 



where ^^^^ = ^(^(t) — ^JJ(t). Hence in our approximation the average radiation reaction force acting on the charged 
point particle a of the system at issue reads 



eaF:i{Xa{T))±l 



h. — 1 ^ 



6=1 



e 2 



_3 

e 2 



ba 



(4 + 2z''{ii + c))e + ie^^z^ ~ 4^"^^^ 



tb I >^"2/ 



(2 - z^ii';, + t,))z, - 2{zPi; + zP's!'^)z, + 2^^ + ^-z% 



\ ^ II. \ -2 Z,, 



3£2 



z ^ -^ z z^ -~\- OZ Z pZ^ 



(21) 



r 



Numerical coefficients at powers of the regularization 
parameter depend on the regularization scheme applied. 
Different powers of the regularization parameter should 
be regarded as independent constants to be taken from 



an experiment. The averaging procedure, which we have 
used, gives only the relations between orders of magni- 
tudes of these constants. 



It is useful to introduce physical masses of particles 

(22) 



eaq -1 

^£ " ■ 
2 



and redefine z(t) by the requirement ((SJ with respect to 
them. Then the equations of motion for the center of 
mass look like 

(M + G)z^ = F^'' + F^6) + F^^(gi- + d") 



-^ per ^ ^ ^fj. 



G = Fp,iPd% 



M 



N 
a=l 



G 



'-^{qQ-d') 



where F^° = f 9^['z^ + z^ip 
and 



(23) 
is the Lorentz-Dirac force 



F,(6) = -r£iq2 



3£* 

8 



(4) 
Z I 



3.. 2.. 

— z z,, 

2 ^ 



■ oz z pZ^ 



(24) 



Formally, the rigid term 1)24(1 is nothing but the counter 
term which, in addition to an ordinary mass renormal- 
ization, compensates the divergencies emerging in six- 
dimensional classical electrodynamics of a point charged 
particle |j, |^ . Such rigid terms have long been known 
for nonrelativistic models of extended charged particles 
in the four-dimensional space-time (see, e.g., 21j and ref- 
erences therein). 

From (|22|) and the condition of positiveness of the total 



^N 



a=i "^a ^^ recover 

q^/2M < e^l"^. 



(25) 



Furthermore, in what follows we imply the severer esti- 
mation 



q^lM < e, 



(26) 



which nevertheless is the reasonable one. 

In order to close the system (|23|l and obtain the effec- 
tive model for a charged object we have to know the time 
evolution of the dipole moment. To this end we make use 
of Eqs. (|14|l . (|21|l again. At lower orders they are 



m„ 



'dr 



^M+C;pit ^eaS^HQ^;~d, 



(27) 



The equations of motion describe the dynamics of me- 
chanical moments therefore we should determine rela- 
tions between electromagnetic and mechanical moments. 



The simplest way to do it is to pick out different species 
of particles. 

We call particles to be of the same species iff they have 
an equal ratio A = e/fii. Let the system be made of K 
species of charged particles with charge to mass ratios A^ , 
s — \, . . . ^K. Then summing in Eqs. H27f) over particles 
of the same species we obtain 



d 



dl + (qs + iJ^dpipJ = XsS 5 (qdl^ - qsdf,) 

+KFpAqsi'' + C) + Xsd1dpF^,z\ 



where qs and d^ are the total charge and dipole moment 
corresponding to the species s. Equations H23() and H28() 
must be solved perturbatively starting from the Lorentz 
equations 



Mip = qFp^z", 



(29) 



and taking into account the estimation (|26|l . In doing so 
we eliminate all the higher derivatives of z^ (r) . 

While the total charge of the object under considera- 
tion is sufficiently small, i.e. q^ /M <^ ^, in substituting 
to (|28|l one can disregard the Lorentz-Dirac force, G-term 
and the rigid term in the equations of motion for the cen- 
ter of mass (|^ and rewrite the equations lP5|l as 

dt + ^Id^^^^-^P-^")^^' = ^^^~Hqd'^, - qsd^) 



gj A, --t ) --^dPFp^z" 



F z" 



+ Fp, ( A.< - ^d^ 



A,d^ - ^d" ) dpFp^z"". 



(30) 



In the case of q^ jM w ^ some additional terms appear. 

Thus we obtain a closed system of equations describing 
in our approximation the evolution of a charged object. 
To put it in another way, singled out essential degrees of 
freedom and averaged over or neglected the rest of ones 
we derive the aforementioned effective model. 

The equations H23|) and H28|) . H3U|) can be regarded as 
the starting point for further perturbation theory in ^ to 
take into account higher multipoles effects. For example, 
these last inevitably appear when one needs to consider 
the self-interaction corrections to H28|) to the same accu- 
racy as in Eqs. H23|l . Consequently, it is worthwhile to 
obtain the time evolution for the second order in ^ mul- 
tipoles. Similarly to the above considerations from Eqs. 
(|27|1 we find at lower orders 



A.T;, + pr^^ d^z,) =\s\e-^ q Q^, - Ql^pr" P^l d^d^ 



^^^(pPp'^) 



'i><,dl--z''[Ql^wl+S;. 



o{p 



F,A„z' 



dT 



Sf. 



%v + Qli^pi'^v) + rf(^^,.) + T.ki^'i' QU p<) +^U) 






= A,U"^ 



2T,! 



dlF^)pi''-^ 



V'^] 



QUp<^+s, 



pu\ 



\s{e '^du4l 



1 






K) + Q^ppa''F,).i- 



F^ 



Qf^p^"^.]-^^ 



where prj;^ 



pr;;, and T^^ = E "^a pr^ pr^^ CpC can (|2Sll, (EHl, (EU for the effective model. 



be interpreted as the intrinsic stress (momentum flux) 
tensor^ corresponding to the species s. As all the intrinsic 
quantities we have defined it is orthogonal to i^. Again 
we arrive at the closed system of equations of motion 

^ Strictly speaking, the stress tensor is JdT(5*(x — z(r))T^i,(r). 
^ This case is of importance since it can be considered as the base 
for perturbation theory in the deviation of the ratios Ai from 



J 



(31) 



Provided that the charged object is composed of parti- 
cles of one species^ the intrinsic dipolc moment vanishes 
due to © and the system of equations (|31|l looks like 

their mean value. 



ATujy — A 



[qe ^ (^Q^^ 



d 

d^ 



'^^lu + W(p 



<3(MPP<i 



AT^fp^p^) 



\^^(q 



p^ pl'(p 



2(p2 



ip<y Pr^.) 



-On 



^pv 



2[p^ 



Qpa Wv] +'5'piy] 



'Sp(f_,] Fj,)„z" 



I (Q(p.pr^ -P"rfpQ.p - ^(Mp) dPF^^^z-], 



A{2ge-iQ^, + 2T^, 



Prf^ Q'yp + Qif^aW^p + 'S'(p 



A r 



Prfp Q'yp + Q[p<TPrp + S[f,p 



F" , 
FP 



Q(pp5''^.).i"}, 
-XQ^^pdPF,],z\ (32) 



r 



The equations of motion for the center of mass (|23|l are 
modified into 

(M + G)zp ^F^"" + F^6) + qF^^z'' + ^QP-dp^F^.z'^ 

+ i(pr'"^ g:; + Q^-p-r;; + sndpF^., 



1 



1 



G =-QP''dpF,,z'^ + -{SP" - prPQ-^Fp., 

(33) 

where we add the first corrections in i^ to the external 
force acting on the charged object as a whole. 

For instance, for the unform external field i^^i, = const 
and constant quadrupole moment Q^y = const the last 
equation of the system (|32|l is nothing but the Bargmann- 



Michel-Telegdi equation 
(|32|l are equivalent to 



MUM- The rest of Eqs. 



AT,,, 



jSi^pK) + W^'Sp,) - Age-iQp.. (34) 



The trace of this equation gives rise to the virial theorem 
in our approximation. Notice as the quadrupole moment 
is constant there is an inertial frame in which Q^^^, be- 
comes diagonal for all times. The condition Qp,yz'^ = 
implies either this is the comoving frame and Zp — or 
at least one of the space-like diagonal elements is zero, 
i.e. it is sufficiently small. If only one diagonal element 
is zero (the case of a thin-plate object) then the vector 
of intrinsic magnetic moment is parallel to the respective 
axis. If two or three elements vanish then S'^j, = 0. 



As we have already noted the terms in Eqs. H21() and 
(|23|l at half-integer powers of e are Lagrangian. In order 
to simplify the Lagrangian we write down it in the natural 
parametrization and keep only the terms which are at 
most linear in the transverse gauge ((BJ: 



N „ _5 



„b=l 



+ ^[(c'")'(i + z^i^p) - 2^'"^;^"^] 



£ 2 

"32" 



[16 + lei^Cp + 8^"^^ + 3(e 



ba\2--2 _ a-pta-aph 



Z 



Si^^i^C 



(35) 

Varying the effective action with the Lagrangian density 
(|S5|l with respect to ^jJ(t) we arrive at Eqs. (|^ . varying 
it with respect to z^^t) and dropping out the terms of 
higher order in ^ we obtaining Eqs. (|23|l . In both cases 
the Lorentz-Dirac force are not of course reproduced. 

To conclude this section we derive the effective equa- 
tions of motion for a neutral system of charged particles 
at the pointlike limit, i.e. at g = and ^ — > 0. Rigor- 



ously, we consider the system for which 

ill < {e^llf, q^f < d\ 



(36) 



where d is a magnitude of the intrinsic dipole moment. 
We also assume the estimation (|18f) is fulfilled for n — I. 

In what follows we are interested in the equations of 
motion for the center of mass of the object in question 
taking into account the first correction due to radiation, 
i.e. to the accuracy of the leading non-Lagrangian contri- 
bution like the Lorentz-Dirac force in the case of charged 
particle. At Fig. ^ the required terms are pictured as 
three dots at the third line corresponding to the order of 
i^. In view of H36|l other contributions are neglible. 

The action for two Lagrangian terms can be evidently 
deduced from H35|) . Its Lagrangian density is equal to 



e 2 



£ 2 



Lg=o = d^+ [4rf^ - 3d^z^ - A{7fdp 
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(37) 



SiPdpS^da 



■ 2{zPdp 



Thus we have to find the non-Lagrangian term only. Te- 
dious calculations analogous to the case of a charged sys- 
tem show up that the non-Lagrangian part of the self- 
force exerting on the system as a whole looks like 



F. 



4 ,2(5) 



p - Y^d z p + -dPdpZp 
2 



(4) 2 r 



2dPdp + (P + 2{zPdpf + i^d^ 



3 1'^'^^ + ^ 



{iPdpf + -i 



2^2 



2 1- • •• ■ •• 

Zp - ■^\d''"dp - {zPdp + iPdpf - 2zPdpZ''da 



+ 2zPdp{-z''d„ - i^'da) - 'i^ W'dp + l{zPdpf 



bzP'Zpd^d^ + d^ 



2 .pi5) ..4 

-z'' z p- z 



'"dp + 'zPdp + ''zPdp + \^fpdp + z^^Pdp - zPdp) 





'z'Pdp - Sz'^zPdp 



dp + ^ 'z^dpdp - -zPdpdp. (38) 



r 



Gathering all the terms together we arrive at the follow- where F^'^s'' denote Lagrangian self-forces coming from 
ing equations for the center of mass: (|37|l . To the same accuracy the time evolution of the 

intrinsic dipole moment obeys 

MZp = Fl^ar _,_ prr ^ p^J. ^ dPdpFp.z", (39) 



dT 



e 2 



d'p + {qs + zPd'p)ip = -Xsqs\ e ''df, + ——[ 2d^ + -z^dp - iPdpi^ - 2iPd_ 



-^P''-p 



iPdpZf, + {zPdp + 2iPdp)zp + {iP'dp + 2z^zPdp)z^ - zP'z pdp - 'i^d^ - d 

+ Xs K.(gsi" + <) + dPdpFp.i'^ 



(40) 



The contributions at the first and third lines are simply 
obtained from Eqs. (|21|l . while the term at the second 
line is derived ab initio. 

So, we derive the equations of motion for the effective 
model for a neutral pointlike object with dipole moment. 
Note that in a similar manner the equations of motion 
can be derived for the effective model for a neutral point- 
like object with magnetic moment and vanishing dipole 
moment. 

It is reasonable to suppose that the energy of the self- 
interaction is much smaller than the rest energy of par- 
ticles in the absence of interaction, i.e. 



dVe3/2 ^ ^.^_ 



(41) 



In that case all the higher derivatives of z^{t) are per- 
turbatively expressed in terms of the first ones by means 
of equations (|39|l . As far as d ^ are concerned the similar 
procedure can be applied to Eqs. (|40|l under the assump- 
tion that XiQiP /e^ are of the order of unity or less. 

Let us extract some physical information from the huge 
equations H39|l . H4U|) . First of all, as one can see from Eqs. 
(|38|l the total radiated power at small accelerations of the 
center of mass of the object is equal to 



' M ^ o dpZfj,, 



(42) 



and coincides with the known nonrelativistic expression 
(see, e.g., [T^) for a dipole radiation. Besides, as it fol- 
lows from Eqs. (|38|l in the nonrelativistic limit^ the total 
radiated power is 



-P,^ 



ld'''dp-^^i2d''lp + dPdp) 



(43) 



Whence the average emitted energy per unit time in the 
lab frame amounts to fdci. 

In the absence of external fields and at zero accelera- 
tion of the center of mass the equations of motion (|S^ . 
(|^ are reduced to® 



d 
d^ 



ae 



IdP'dp, 



-as ^dn 



2a- 



(44) 



respectively. In order to the system of equations 144(1 
possesses nontrivial solutions with the characteristic scale 



'^ Recall that on recovering the velocity of light every overdot con- 
tains 1/c. 

* The second equation in 1441 resembles the equation of motion for 
a damped linear oscillator. The existence of self-oscillations of 
charged distributions has been known long ago (see, e.g., l24il ) . 



of variations more than or equal to / the positive quantity 
a — J2s=i ^s<ls is supposed to meet the condition 



a'^/e <£2/; ^ i. 



(45) 



The system (I44|l can be satisfied by the second equation 
only provided the additional requirement 



const 



(46) 



is fulfilled. 

By the general procedure we must solve Eqs. (|44() per- 
turbatively, but in our case we can do it exactly and make 
approximations directly in a general solution. The char- 
acteristic numbers associated to the second equation in 
(|44|l are rather complicated and we write down them in 
an approximate form 



Ai,: 



±iuj 



1 ■ 






27' 



2a 



> 2 + ae i 

(47) 
taking into account the estimation (|45|) . The solu- 
tion corresponding to the third characteristic number is 
known as runaway solution. We should drop it out since 
it does not meet the initial hypothesis concerning the 
characteristic scale of variations of the fields Cu (''')) i-^- 
A3 ^ l"^. By the assumption (|45|l the damping factor 7 
is small, 7 <IC /~^, that is the reason that the solution of 
the second equation in the system l(^ 

^^(t) = [a^ cos (cjt) -t- 6^ sin (wT)]e"''''^, 

2_,2 .n r. (48) 



^h' 



UpbP = 0, 



approximately satisfies the requirement H46|l . This solu- 
tion represents a freely moving slowly rotating dipole. 



IV. HYDRODYNAMICAL APPROACH 

In this section we regard a Poincare-invariant hydro- 
dynamical approach to the description of a system of 
charged particles thought as a relativistic perfect fluid. 
We briefly reformulate basic definitions previously intro- 
duced within the particle framework and prove the equiv- 
alence of the hydrodynamical and particle approaches to 
the problem of radiation reaction for multipoles under 
some assumptions given below. In conclusion of this sec- 
tion we give a generalization of the notion of multipole 
moments to extended relativistic objects (branes). 

A simple general covariant generalization of the ac- 
tion ^ to the hydrodynamical case can be constructed 
as follows. Suppose given a 3-brane N with coordinates 
{r*}, i = 0, 1, 2, 3, that is embedded by a diffeomorphism 
x^{t) into Minkowski space R'^'^. If the brane N is the 
space-time itself one can think about such a diffeomor- 
phism as a general coordinate transformation in it. Let 



us introduce p^{t) and e*(r) that are vector densities on 
the brane N and describe a mass flow and an electric 
current respectively. Then an obvious generalization of 
the action (Q takes the form of^ 



S[x{t),A{x)]=- / d\Jp-fPh, 

J N * 



<fx 



^^'' + 16^^'-^'' 



(49) 



where hij — diX^djX^ri^^y is the induced metric on the 
brane iV, it being flat inasmuch as the space-time is flat. 
The vector densities p^ir) and e*(T) are supposed to van- 
ish at spatial infinity. The electric current density on the 
target space (the space-time) is 

f(x)^ [ d'^T6\x - x{T))e'{T)d^x^'iT). (50) 

Jn 

Of course, the (5-function can be integrated out resulting 
in the Jacobian, but we leave it intact to keep an analogy 
with |(2Jl. In other words the matter action in (|49|l defines 
the dynamics on the infinite-dimensional group of dif- 
feomorphisms^" x{T'^,f) of the space-like hypersurfaces: 
T^ = const and its image in the space-time. Hereinafter 
we denote f — {t^,t^,t^}. 

An invariance of the action functional H49|l under gauge 
transformations of the electromagnetic potentials implies 
the charge conservation law 



d,e' = 0. 



(51) 



Then the equations of motion for the matter obtained 
from the action (|49|l read 



5,: 



P'p'djX^ 



j-^M 



-T^jyC OiX . 



The matter energy-momentum tensor is^^ 

TZt - Id'r6\x x{T))^d,x^d,x\ 



(52) 



(53) 



^ The action functional for th e co upling of branes to (non)-Abelian 
gauge fields is suggested in l25l . 

^^ A group-theore tica l approach to hydrodynamics can be found, 
for example, in l26l , where the equations of motion of a fluid are 
regarded as the equations of geodesies on the infinite-dimensional 
group of diffeomorphisms generated by the fluid flows, which is 
equipped with a right-invariant Riemannian metric. 

^^ Notice that the matter action I49i can be easily generalized to the 
case of a perfec t fluid with pressure (see for other Lagrangians 
|271 l28l 1291 l30l and references therein). Introducing a scalar 
density p(t) on the brane A'^ we define the action functional as 



S{x(r)\ 



\J p^ -I- p In 



-h 



where h is the determinant of the induced metric, the vector 
density P^(t^ is the internal energy flow and "pir) characterizes 
the fluid pressure. 



that is the energy- momentum tensor for a relativistic 
dust (see, e.g., |31|). 

In a certain sense the hydrodynamical approach is re- 
duced to the particle one under the assumption 



e\r)^\{T)p\r), 



(54) 



where A(t) is some scalar function on A'^ describing the 
charge to mass ratio. If that is the case Eqs. (|52|l give 
rise to 



d^p' = 0, p'9,A = 0, 



(55) 



i.e. to the mass conservation law. Provided the require- 
ment ()54(l is fulfilled the equations of motion (|52(l possess 
a "partial" reparametrization invariance, which implies 
orthogonality of the equations to p^dix^ . 

Now bearing in mind the assumption (|54l) we proceed 
to multipoles. Let z^(jj) be a naturally parametrized 
worldline in Minkowski space, then similarly to lO we 
claim that 



ip{r^)i''{r) = 0, 



i'(r°) - 1, 



(56) 



where C^(t) = x*^[t) — z^ir^) and overdots denote the 
derivative with respect to r". The condition (|56|) par- 
tially specifies a parametrization on the brane N and 
foliates it into a family of the space-like hypersurfaces 
T^ = const. Thereby the condition (|56|l fixes the gauge 
and spoils the reparametrization invariance of the equa- 
tions of motion. The trajectory z^(t°) of the center of 
mass is defined by 



dfp°{T)C{T)=^0. 
The definition of multipoles © transforms into 

rfre°(T)^^,(T)...e^,.(T), 



(57) 



dfe\T)^^Ar)■■■^[^.„_Ar)^^Ur)w'^^y 



(58) 



It is easy to see that the Lienard-Wiechert potentials 
are expressed in terms of the multipole moments (|58|) and 
their derivatives. Indeed, expanding in powers of £, the 
(^-function in the formula ij^l adjusted to the hydrody- 
namical case we obtain expressions containing ^^(t) in 
the following way 



df,5('=)((a;-z(r"))2) 
X ^,dr) ■ ■ ■ eM„-i(^)e'(T)a,(z,„(T°) + e^Jr)), 



(59) 



that are obviously expressed in terms of the multipole 
moments (|58|) . 

As to mechanical moments are concerned the equations 
of motion for particles ((SJ suggest a general correspond- 
ing rule between mechanical moments in the particle ap- 
proach and hydrodynamical one 

^ ^ /df, m„ ^ p°, d/dr ^ p'/p^di. (60) 
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In our case the ambiguity of this rule must be settled 
by the prescription that, at first, masses nia are removed 
from under the differentiations with respect to r. For in- 
stance, the intrinsic stress tensor arising in H31|l is rewrit- 
ten as 



r„. 



-PR 



dr—^di^pdj^a pr^ pr;^. 



(61) 



In the sequel we do not repeat calculations of the pre- 
ceding section, but we show that the resulting equations 
of motion of the effective model are identical to ones de- 
rived in the particle framework. Besides, in spite of this 
assertion the equations of motion for the center of mass 
(|23|l were rederived by direct calculations. 

Similarly to the considerations regarding the Licnard- 
Wiechert potentials it can be proved that, firstly, the 
equations of motion for the center of mass are expressed 
in terms of the mechanical moments and multipoles (|58|l : 
secondly, on introducing species^^ the evolution equa- 
tions for the intrinsic species moments both mechani- 
cal and electric (magnetic) are expressed in terms of 
the intrinsic species moments. In the hydrodynamical 
framework the species are domains on the hypersurfaces 
T° = const of the constant ratio A(T°,f). The form of 
these domains depends on t" to be consistent with (|55|l . 

Whereas the equations of motion and expressions for 
multipoles in the hydrodynamical case pass into respec- 
tive expressions for the particle case with the assumption 



p^(r) = (p"(f),0). 



(62) 



we infer that the above-mentioned evolution equations 
are identical to the equations derived in the preceding 
section. Notice that the mass renormalization (I22|l in 
the hydrodynamical framework looks like 



P 



e'g 



^^p' = 0. 



(63) 



In conclusion we point out that the elaborated ap- 
proach to the definition of the multipole moments is a 
quite general one and can be used to describe the mul- 
tipole moments for a charged fluid or system of charged 
particles approximated by not only a point particle but 
also a brane, i.e. to define in a Poincare-invariant way 
the intrinsic linear density of dipole moment for a string 
or the intrinsic surface density of magnetic moment for a 
membrane etc. 

Let us split the coordinates {r'} on the brane N into 
{r°,f°}, T° is included to {r°}. Suppose given a brane 



^^ We emphasize once more that separation into species is the sim- 
plest means to obtain a system of ordinary diflFerential equations 
from the integro-differential one reducing to the multipole mo- 
ments I58i only. For example, another interesting choice of the 
ratio A is to assume it to be the kernel of a differential operator 
compatible with 1551 . 



M with coordinates {cr°} that is embedded into the 
space-time by a smooth mapping z^{a). The induced 
metric on the brane M we denote by hat- Suppose the 
brane N is also submersed into M by the mapping 



ip: N ^ M, 



ip-[T,T) 



(64) 



i.e. the mapping ip matches systems of coordinates on 
the branes TV and M. For our choice it merely identifies 
the coordinates {t"'} and {(7°'}- On the one hand this 
mapping foliates the brane N into a family of the space- 
like surfaces r" = const. On the other hand it foliates the 
brane N and, consequently, the space-time into a family 
of surfaces diffeomorphic to the brane M. Then we claim 
that the vector field {p^p'^hi,c)~^ p°'da is projectable onto 
leaves of the latter foliation, i.e. 



[d-a,^d,]=0, (65) 



hereinafter p^(t) — p^{T)p''{T)ht,c{a). Evidently this re- 
quirement can be satisfied by appropriately chosen p* (t) 
and x^{t) for any given vector density 

V^'ix)^ j d'^T5^{x-x{T))p'{T)d,x^'{T), (66) 

JN 

on the space-time. 

The transverse condition H56|l becomes 

daZp{cj)iP{T) = 0, (67) 

where ^^(r) — x^{a, f) — z^{a). This condition partially 
fixes a system of coordinates on the brane N, i.e. it 
specifies the space-like surfaces t° = const on the brane. 
The center of mass condition modifies into 



df^^{^e{T)=0 ^ dfp^{T)^^{T)^0, (68) 



and the intrinsic electric and magnetic multipole mo- 
ments are defined by 



dfe^{T)^^Ar).--it.Jr), 



rf^eXr)e^,(r)...C[^_,(r)5,ep(r)pr^^j, 



(69) 



where pr ^ — S'j^ — p'^p^daz'^dbz^j^/ p^ and we assume the 
relation (|54|l is fulfilled. The mechanical moments are 
generalized in an obvious manner. Clearly, the so defined 
multipoles both mechanical and electric (magnetic) are 
orthogonal to daZ^ and zeroth multipoles are conserved 



daJdfe''{T)^daJdfp-{T)=0, 



as long as corresponding currents are conserved. 



(70) 
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V. CONCLUDING REMARKS 

Let us summarize the main results of our research. We 
have investigated the effective dynamics of a system of 
charged particles within two approaches - the particle 
itself and hydrodynamical one. We have described the 
effective dynamics of such an object by means of a sys- 
tem of ordinary differential evolution equations for the 
intrinsic multipole moments and the center of mass. We 
have proved the equivalence of two examined approaches 
to the problem of radiation reaction for multipole mo- 
ments. In passing we have derived the effective model 
for a neutral pointlike system of charged particles. 

These results can be extended to several directions. 
We point out some of them only. It would be interesting 
to study the effective dynamics of a string or membrane 
with intrinsic higher multipole moments like that is given 
in [l2| for an electrically charged string. Another direc- 
tion for further research could be in a study of a general- 
ization of the effective dynamics of a high-current beam 



of charged particles to include the matter term proposed 
in Section IIVI Other points could be the nonrelativis- 
tic effective dynamics of a neutral pointlike object with 
dipole moment or may be the effective dynamics of a neu- 
tral pointlike object with magnetic moment and vanish- 
ing electric dipole moment subjected to some simplifying 
assumptions. 
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